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ABSTRACT 


(wn: 


Percentage points of a new distribution involving a confluent- 
hypergeometric distribution obtained by Khatri and Rao (1985)-~are 
tabulated. The use of the tabulated values in obtaining a lower 


confidence bourd for the realized signal to noise ratio based on 
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rs an estimated aiscriminant function for signal detection is ex- j i. 
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ale 1. INTRODUCTION 

Las 

ete Let X be a random variable having one of two possible p-vari- 

ae ate normal distributions N (0,2) and W (6,2) when X is real (or 

i N (0,2) and ane when X is complex). In the sequel, we follow 

ey 

ia 

ue the practice of giving the results for the real case first, and 

Pa | 
ty then for the complex case within brackets. For exact expressions | 
: of the p-variate normal density in the real and complex cases ref- 

soe erence may be made to Khatri and Rao (1985). We consider the case 

at) = 

eo where 6 is known and £ is unknown but an estimate f 1s of £ 
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é | : 
et is available, where S has the Wishart distribution WCE) (or W(£.2)) 
of on £ degrees of freedom. In such a case, the estimated linear dis- 
ky, 

s criminant function is proportional to sts71x (or ses), The 

ct discriminatory power or the probability of correct classification 
Bt j when the estimated discriminant function is used on future samples 
m is a monotonic function of the signal to noise ratio 

AA es ee =1 

a4 p(s,t) = SO | or pis,z) = SS) ay) 
aA 6's “£8 6 6é*§ “rs “6 

“ which involves the unknown £. By the Cauchy-Schwartz inequality, 
- (1.1) is < 6tztg (or sepa) , the signal to noise ratio when £ 
ot 
LN is known so that there is some loss of information in using an esti- 
& 

mated discriminant function. The problem is to make an inferential 

o: 

ay statement on the realized ratio (1.1) in terms of known quantities 
4 

: 6, S and f. This can be done in several ways as shown in the next 
oe 
eo. section. 

ts 2. INFERENCE ON SIGNAL TO NOISE RATIO 

“e ened 

+ 

5 The key result for this purpose is Theorem 1 in Khatri and Rao 
. '(1985) where it is shown that in the real case 
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a (6's7 78)? girate (2.1) 
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1% In the complex case, defining B and G with 6' replaced by 6* 
7 in (2.1), it is shown that B and G are independently distributed, 
x 
4 with the p.d.f. of B as 
a T(£+1) p(f-pHL) (5 _5)p-2 (2.4) 
I (£-p+2)T (p-1) : 
i ~ 
ng and that of G as 
"4 
L) 
N 1 f 
fone ek BB: 7 EP 
t fpr? 8 (2.5) 
fl The distribution (2.4) was earlier obtained by Reed, Mallet and 
Y Brennan (1974), The distributions (2.3)-(2.5) are independent of 
‘ the unknown parameters which enables inferences to be drawn on 
t (1.1) through the pivotal statistics B and G (or B and G). 
v 
‘s 1) Using the expressions for the moments of the beta distribution 
‘a (Rao, 1973, p. 168) 
a 
= f-pte 

9 E{B(or B)] 41 
iS 
t 
r which gives 
Ks : 
4 E[p(S,f)(or p(S,z))] = ee ("27 "sor 6*2718)}. (2.6) 
., 
If the average efficiency is to be maintained at about half the 
¥ 
‘s optimal efficiency, then from (2.6) we have 
3 $optt oboe pa 
a f41 5 or £ = 2p (2.7) 
3] . 
! for both the real and complex cases. The result for the complex 
2 case is mentioned in Reed, Mallet and Brennan (1974). Similarly 

; we can equate the ratio in (2.7) to any desired ratio other than 
i (33) and obtain an expression for the degrees of freedom f needed 
‘ for the estimation of Cf. 
* 2) Perhaps a more satisfactory way of using the distributions (2.2) 
2 and (2.4) is as follows. Let b (or b.) be the lower a% point of 
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oy the distribution (2.2) (or (2.4)). Then we can make the confidence 


statement that 


p($,Z) > b 6'Z 8 (or o(S,2) > b,6*E7 6) (2.8) 


Ny —_ with a confidence coefficient of (1-a)%. : 

3) The results (2.7) and (2.8) still involve the unknown quantity 
=. We raise the question whether the actual magnitude of o(S,Z) 
for given S can be assessed through known quantities. Using the 
joint distribution of B and G as in (2.2) and (2.3), or B and G 

as in (2.4) and (2.5), it is shown in Khatri and Rao (1985) that 
the estimate 


i a _ (£-pt2) (f-p-3)_ 2 a _ (é-pt2) (f-p-1) 
f 0 (8,2) €(£+1) D (or 6(8,2) £(£+1) DE *) 
ie 2 1 1 Min 2,2 
acy where D. = £5'S & (or £6*S 6) has the property - E{e(S,2) = gDo1” 
aS 
‘ey 
‘See 
2 a 
=F [0(S,2) = p(S,2))* (or = E [B(S,2) - 8(S,2)]%). 

ay 

pe. a a 
Nt: so that p(S,£) (or B(S,£)) is a predictor of 9(S,Z) (or §(S,2)) 

ay 4) Khatri and Rao (1985) also obtained an exact confidence bound 
v for 9(S,£) for the computation of which we provide extensive tables 
ret in this paper. 
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% We define the random variable 
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which has the confluent hypergeometric distribution with the proba- 
bility density function 
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ae as described in Khatri and Rao (1985, Theorem 3 and Remark 5). If 
wt zy (or z.) is the lower a% point of this distribution, then 

at 

a 

Xt ‘ 

i) 

Cy 2z,, 2 i zy 2 

% p(S,2) ar a (or p(S,2) oF Dd) 

“3 

sy provides the lower confidence bound to the signal to noise ratio 
Y p(S,Z) (or p(S,£)) with a confidence coefficient of (l-a)%. Tables 
‘ > 

“ 1-5 given in the next section give the values of 22 (or 2) for 


a. various combinations of p and f, and a = 0.05, 0.25, 0.50, 0.75 
Po and 0.95. 
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i 3. PERCENTAGE POINTS AND SOME APPROXIMATIONS 

f 

Tables 1-5 give the lower a% (for a% = 5, 25, 50, 75 and 95) 
a points of the distributions of 2z and z for different-values of 
At 
ny p and f. Actually zy is obtained as a solution to the equation 
4 
NS 7 
1 (£41) z/y ; 
r 2 (f-p) /2 (p-3) /2 1 -x _(f-p-1) /2 

a = ca VEEN IFES, ERS | (l-y) dy = x dx 
wi : 
iS 
5) 4 
‘ng and z_ to the equation 
: = z/y 
", T(£+1) f-p+1 p-2 
Mn as\ —— iy (1-y)P~“ay 1 -x _f-p 
a} is T(p-1) PF (£-p+2) 0 T(E-ptl) ex ‘dx. 
an 
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‘ To obtain (1l-a)% lower ne cea we multiply the Mahalano- 
re bis distance DF = £5'S” 1s (or £5*S_ Es) by 22, /£ (or z q/f)- 
A 
on We give several approximations to the distributions of z (or 
Ss: z) from which fairly approximate values of ze (or z,) can be easily 
‘ 
“ obtained if p/f is not large. 
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(i) Gamma approximation 
(a) We consider the statistic 


g(p,f)z (or g(p,£)2) (3.1) 


and approximate its distribution by a gamma atekesiut ion G(1,v) 


or G(1, v)). 
To determine g(p,f) and v we equate the first two moments of 
g(p,f)z with those of G(1,v). The equations are 


E[ g(p,f)z] = g(p,f)E(z) = v 
v(s(p,£)2] = (g(p,£)]2v(z) = v 
which give. 


E(z) (£+1) (£+3) 
g(p.£) = 702) ™ (Eni) (£43)-p(pel) 


xs LT E(2) J? _ f-ptl __ (£+3) (£-pt2) ; 
V(z) 2 = (£+1) (£+3) -p(p-1) 


Similarly 
-(p.#) = EZ) 2 (f+) (£42) 
Bee ya) CEPT) #2) -p(p-D) 
~ _ _[E(z)]* _ (£42) (£-p+l) (£-p42) 
eG) (e+) C42) -p(p-1) 


(b) We consider the statistic 
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g(p,f)2° (or g(p,f)2° (3.2) 


and approximate its distribution by a gamma distribution G(1,v) 

(or G(1,v))- To determine g(p,f), ¢ and v, we equate the first three 
moments of g(p,f)z° with those of G(l,v). Taking m = (f-ptl)/2, 

a = (f-pt2)/2 and b = (p-1)/2, the equations are 
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v = g(p,£)T (mtc) TP (atc) T (atb) /{P (mT (a)T (atb+c)}, 

vtl = g(p,f)P (m+2c)l (at2c)l (atbtc) /{T (mtc) T (ate)? (atb+2c) } 

: v+2 = g(p,f) I (mt3c) IP (at3c) IF (atbt2c) /{T (m+2c)T (at2c) I (atb+3c) }. 


These equations give 


4 
x T (m+3c) 0 (at3c) Pf (atb+2c) T (mtc) T (atc) I (atb) | 
i. T(mt2c)P(at2c)P(atb+3c) © T(m)T(a)T (Catbt+c) 
bh 
=2 I (m+2c)T (at2c)T (atb2c) 
y T (mtc) I (ate) T (atb+2c) 
ie from which a solution for c is computed using an appropriate com- 
7 
b> puter program. 
c Similarly, for the complex situation, we take m = f-p+tl, 
Q a = f-p+2 and b = p-1, and determine the values of c, g(p,f) and v 
{ as above. The values of these constants for the real and complex 
in cases are given in Table 6 for some values of p and f. 
nt We find by comparing with actual values that the approximation 
” given by (3.2) is more accurate than that given by (3.1), and the 
b approximation is fairly accurate even for small values of f and p. 
a (ii) Normal approximation 
y (a) Using the Wilson-Hilferty's approximation as modified 
u by Konishi (1981), we have 
Pal 
ww (ov) 8 (BB@22)1/3 _ 14 Sy - NC0,1) 
2vu 9v 
% Lor (9v) 2 (B24) 1/3 yy Ly - n0,1)1. 
a8 2v ov 
reed 
= where g(p,f) and v (or g(p,f) and v) are defined in (3.1). 
x (b) As in (a), we can take 
s : 
re (oy) '#¢(Be22 1/3 4 hy ~ won) 
Ks 2v 9v 


ms i ~¢c 
{ or (9y) 4 (Bieafz-) 1/3 -1l+ a ~ N(O,1)] 
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eae where g(p,f), v and c (or g(p,f), v and c) are defined in (3.2). 
eal The normal approximation is fairly accurate even for small values 
te . of f and p. 
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